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ABSTRACT 

We use high-resolution spectral emission line data obtained by the SERTS instrument during three rocket flights to 
demonstrate a new approach for constraining electron densities of solar active region plasma. We apply differential 
emission measure (DEM) forward-fitting techniques to characterize the multithermal solar plasma producing the 
observed EUV spectra, with constraints on the high-temperature plasma from the Yohkoh Soft X-ray Telescope. 
In this iterative process, we compare line intensities predicted by an input source distribution to observed line 
intensities for multiple iron ion species, and search a broad range of densities to optimize y 2 simultaneously for the 
many available density- sensitive lines. This produces a density weighted by the DEM, which appears to be useful 
for characterizing the bulk of the emitting plasma over a significant range of temperature. This “DEM-weighted 
density” technique is complementary to the use of density-sensitive line ratios and less affected by uncertainties 
in atomic data and ionization fraction for any specific line. Once the DEM shape and the DEM-weighted density 
have been established from the iron lines, the relative elemental abundances can be determined for other lines in 
the spectrum. We have also identified spectral lines in the SERTS wavelength range that may be problematic. 

Key words: Sun: abundances - Sun: corona - Sun: fundamental parameters - Sun: UV radiation - Sun: X-rays, 
gamma rays 
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1. INTRODUCTION 

Plasma temperatures, densities, and elemental abundances 
are crucial parameters required to understand and model the 
solar corona. Despite the fact that differential emission mea- 
sure (DEM) modeling has been available since the 1960s, most 
temperature measurements have been made under the assump- 
tion that the emitting plasma was isothermal. Recent results, 
however, imply that this simplification may not apply to as 
many coronal structures as we once thought. Densities, n e , are 
usually found from the ratio of intensities of two spectral lines 
of the same ionization stage of a given element. When an ap- 
propriate pair of lines is not available, however, densities are 
sometimes estimated from observations of the emission mea- 
sure, EM = n 2 e dl, where the emitting length (/) is determined 
from an inspection of the image and an assumption about the 
filling factor. Elemental abundances affect cooling timescales, 
energetics, pressure balance, and also determine the shape of 
the radiative loss function. They are tracers of the origins of the 
solar system, and have been used to track various components 
of the solar wind, solar energetic particles (SEPs), and coronal 
mass ejections back to the solar source. 

Spectroscopic (Schmelz et al. 2012), SEP (Reames 1995), 
and solar wind (von Steiger et al. 2000, 2010) data suggest 
that the coronal-to-photospheric abundance ratios of elements 
with low first ionization potential (FIP < 10 eV) seem to 
be enhanced relative to those with high FIP (> 1 1 eV). This 
fractionation probably results from a separation of ions and 
neutrals. Two different empirical models have been suggested 
in the literature: (1) low-FIP elements may be enhanced by 
about a factor of four with respect to their photospheric values 
while high-FIP elements are the same in the corona and the 
photosphere (Feldman 1992), or (2) low-FIP elements may be 


the same in the corona and the photosphere while high-FIP 
elements are depleted by about a factor of four with respect 
to their photospheric values (Meyer 1985). However, many of 
the elemental abundance measurements made in the 1980s and 
1990s and used to support one or the other of these models were 
derived with old atomic data and under the assumption that the 
coronal plasma was isothermal. 

In this paper, we take advantage of the improved understand- 
ing of the thermal structure of the corona, better multithermal 
analysis techniques, and most importantly, updated atomic mod- 
eling and ionization balance calculations for coronal emission 
lines. The spectral data from three different active regions ob- 
served during the 1993, 1995, and 1997 flights of the Goddard 
Solar EUV Rocket Telescope and Spectrograph (SERTS-93, 
SERTS-95, and SERTS-97) were described in detail by Brosius 
et al. (1996, 1998, 2000). We have used these spectra along 
with updated atomic physics data from CHIANTI v.7 (Dere 
et al. 1997; Landi et al. 2012) to derive the DEM for the ob- 
served active region plasma. Broadband X-ray filter data from 
the Yohkoh Soft X-ray Telescope (SXT) were used to constrain 
the high-temperature end of these DEM distributions. These 
results were used to derive a DEM-weighted plasma density 
and to test the new set of “recommended coronal” elemental 
abundances proposed by Schmelz et al. (2012), an update of the 
“hybrid” abundances of Fludra & Schmelz (1999). 

2. OBSERVATIONS 

The SERTS (Neupert et al. 1992) payloads were launched 
on Terrier-boosted Black Brant rockets from White Sands, 
New Mexico. Each incorporated a grazing incidence Wolter 2 
telescope, which formed a real image of the Sun on the entrance 
aperture of a quasi-stigmatic spectrograph. The aperture allowed 
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spectra to be obtained along a long narrow slit which connected 
two rectangular lobes, within which spectroheliograms were 
imaged. In our analysis, SERTS images were used only for co- 
registration of the SERTS and SXT data; otherwise, only the 
SERTS slit spectra of the active region targets were considered. 

In each of the three flights (see Brosius et al. 1996, 1998, 2000 
for details), the instrument included a multilayer-coated toroidal 
diffraction grating which enhanced the instrument sensitivity 
over a portion of the spectrum as compared with that of a 
standard gold grating (such as that used in the 1989 SERTS 
flight). SERTS-93 and SERTS-97 used the same multilayer- 
coated grating, which enhanced the first-order response around 
300 A by up to a factor of nine over a gold grating. SERTS-95 
used a different grating, which enhanced the sensitivity of 
the second-order wave band (171-225 A), producing greatest 
sensitivity in the range ~ 182-208 A. 

SERTS-93 was launched on 1993 August 17 at 1800 UT 
and recorded spectrographic data between 100 and 492 s after 
launch. Spectra were obtained from 235 to 450 A; the measured 
instrumental FWHM resolution varied in a known way along the 
slit and was typically ~55 mA; spatial resolution was ~5 arcsec. 
An active region spectrum with 65 prominent lines was obtained 
from an average of the 163 arcsec portion of the slit that bisected 
Active Region (AR) 7563 (at heliographic coordinates SOI 
W15); the slit was roughly parallel to the direction of the loops. 

SERTS-95 was launched on 1995 May 15 at 1800 UT and 
recorded spectrographic data for nearly 7 minutes. The FWHM 
spectral resolution was ~30 mA in second order (171-225 A) 
and ~55 mA in first order (235-335 A); the spatial resolution 
was ~5 arcsec. A spectrum was obtained for AR 7870 (at 
N09 W22), by averaging the spectra over a 226 arcsec portion 
of the slit that bisected the region; in this case, the slit was 
approximately perpendicular to the direction of the loops and 
crossed near the loop tops. Over 100 lines were identified in the 
spectrum. 

SERTS-97 was launched on 1997 November 18 at 1935 UT 
and recorded spectrographic data for nearly 400 s. The FWHM 
spectral resolution was 115 mA; the spatial resolution was 
~7 arcsec. A spectrum was obtained for AR 8108 (at N21 E18) 
by averaging over a 221 arcsec segment of the slit covering the 
region. The spectrum included 95 emission lines — 52 identified 
spectral lines plus 43 weaker features that correspond to blends 
or have uncertain identification. 

The SXT images are taken in broadband soft X-rays 
(approximately 2-60 A). Details of the instrument and opera- 
tion modes can be obtained from the paper by Tsuneta et al. 
(1991); the Yohkoh Analysis Guide (www.lmsal.com/SXT/) 
documents much of the available software for data reduc- 
tion and analysis. Also, please see the Yohkoh Legacy data 
Archive (ylstone.physics.montana.edu/ylegacy/). SXT provided 
cotemporal high-temperature (>2 MK) images of the targeted 
active regions. The data were co-registered by the SERTS 
team. For SERTS-93, they used partial-frame, high-resolution 
(2.455 arcsec pixel) images in three filters: thin Al (379.95 
DN s" 1 ), AlMgMn ( 1 92.8 1 DN s~ 1 ), and thick Al (9.00 DN s“ 1 ). 
For SERTS-95, they used full-frame images with 4.91 arcsec 
pixels in two filters: thin Al (2968.10 DNs -1 ) and AlMgMn 
(15.73.45 DNs" 1 ). And for SERTS-97, they used full-frame 
images in the AlMgMn filter (821.57 DN s -1 )- 

3. ANALYSIS 

In this paper, we analyze active region data obtained simulta- 
neously by SERTS and SXT. The SERTS iron lines used in this 


analysis are listed in Table 1 and the non-iron lines are listed in 
Table 2. These data, including the line identifications and wave- 
lengths, are from the original SERTS papers by Brosius et al. 
(1996, 1998, 2000). If a line appeared in more than one data 
set, we listed the wavelength from Brosius et al. (2000). The 
uncertainty of the line centroid is typically less than 10 mA for 
most lines (Brosius et al. 1996). We have also listed the wave- 
lengths from CHIANTI v.7 for comparison. By themselves, the 
available SERTS data do not constrain higher temperature ac- 
tive region plasma. The SERTS wavelength range does include 
higher temperature flare lines, but none of them were detected 
for these active regions, and the signal-to-noise level was not 
high enough to use the absence of the lines as a significant high- 
temperature constraint. Combining SXT data with the SERTS 
data provides the constraint on the high-temperature end of the 
DEM curve which is typically poorly defined in studies using 
EUV data alone. 

Although the majority of the SERTS emission lines have sig- 
nificantly cooler temperature response than that of the dominant 
SXT response, there is sufficient overlap of temperature sensi- 
tivity to permit consistency checks between the SERTS and SXT 
analyses. So, the combined data set is able to yield a reliable 
multithermal analysis over the range (1-10 MK) needed to ade- 
quately characterize the thermal distribution of coronal plasma 
in the observed active regions. 

The SXT response functions are computed by folding a 
synthetic solar spectrum through the measured effective area 
for each of the five different filters. Standard SXT analysis 
was done with the ionization balance calculations of Arnaud 
& Rothenflug (1985) and Meyer (1985) coronal abundances. 
Standard SERTS analysis, however, used Arnaud & Raymond 
(1992) ionic fractions for iron and Feldman (1992) abundances. 

Schmelz et al. (1999) showed that the choice of abundance 
normalization and/or ionization fraction used in the SXT analy- 
sis can severely affect models of the coronal plasma that require 
the temperature or emission measure distribution. Given the sen- 
sitivity of the SXT filter responses to the change in elemental 
abundance normalization and iron ionization fractions, joint use 
of the SXT and SERTS data requires that the SXT response 
functions and the SERTS contribution functions be generated 
with the same set of elemental abundances and ionization frac- 
tions. For the multithermal analysis presented in this paper, we 
adopted the ionization fractions of Bryans et al. (2009), which 
the CHIANTI team has recalculated with a temperature resolu- 
tion of log T = 0.05. For our initial best guess of the elemental 
abundances, we have adopted the recommended coronal 
abundance set of Schmelz et al. (2012). 

3.1. Method 

We use two different DEM methods in our multithermal 
analysis. The first, DEM_manual (Schmelz et al. 2011a, 201 lb), 
is a forward-fitting technique with a manual manipulation of the 
DEM. The best fit is determined from a / 2 minimization of 
the differences between the observed and predicted line fluxes. 
The main advantages of this method are that no smoothing is 
required beyond that imposed by the temperature resolution 
of the spectral line contribution functions or SXT response 
functions. Also, no a priori shape (Gaussian or double Gaussian, 
for example) is imposed on the final DEM curve. The main 
disadvantage of DEM_manual is that it is not usually possible 
to explore a broad parameter space, so families of solutions 
might be missed. 
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Table 1 

SERTS Iron Line Intensities (erg cm~ 2 s -1 sr 1 ) 



Ion 

log T 

SERTS 
A. (A) 

CHIANTI 

A (A) 

SERTS-93 

SERTS-95 

SERTS-97 

Notes 

1 

Fe vm 

5.70 

185.221 

185.2133 


305 ± 42.3 


a 

2 

Fe vm 

5.70 

187.247 

187.2407 


15.6 ± 7.40 


b 

3 

Feix 

5.90 

171.074 

171.0733 


208 ± 62.0 



4 

Feix 

5.90 

218.943 

218.9373 


98.9 ± 30.1 


c 

5 

Fex 

6.05 

174.526 

174.5310 


654 ± 83.0 



6 

Fex 

6.05 

175.265 

175.2630 


135 ± 32.4 



7 

Fex 

6.05 

177.240 

177.2400 


335 ± 47.5 



8 

Fex 

6.05 

184.534 

184.5370 


164 ± 23.8 



9 

Fex 

6.05 

190.046 

190.0370 


47.9 ± 9.40 



10 

Fex 

6.05 

345.735 

345.7380 

106 ± 12.7 


57.0 ± 9.60 


11 

Fexi 

6.15 

180.382 

180.4012 


3440 ± 341 


a 

12 

Fexi 

6.15 

181.125 

181.1302 


123 ± 35.6 



13 

Fexi 

6.15 

182.166 

182.1672 


237 ± 33.4 



14 

Fexi 

6.15 

184.800 

184.7930 


31.5 ± 14.0 


b 

15 

Fexi 

6.15 

188.214 

188.2165 


385 ± 77.3 



16 

Fexi 

6.15 

188.295 

188.2994 


276 ±51.6 



17 

Fexi 

6.15 

189.185 

189.1234 


32.8 ± 8.90 


c 

18 

Fexi 

6.15 

189.722 

189.7112 


26.4 ± 6.70 



19 

Fexi 

6.15 

192.813 

192.8137 


68.3 ± 9.80 



20 

Fexi 

6.15 

202.415 

202.4242 


73.7 ± 17.3 


c 

21 

Fexi 

6.10 

308.548 

308.5445 

28.9 ± 11.4 

120 ± 37.4 

33.8 ±9.6 

b*,c* 

22 

Fexi 

6.15 

308.991 

308.9916 



29.1 ± 11.3 

b 

23 

Fexi 

6.10 

341.129 

341.1136 

71.1 ±9.23 


37.9 ± 6.60 

a 

24 

Fexi 

6.10 

349.032 

349.0466 



7.5 ± 2.50 


25 

Fexi 

6.10 

352.674 

352.6709 

177 ± 20.2 


136 ±22.0 


26 

Fexi 

6.10 

356.553 

356.5193 

19.3 ± 8.38 



b 

27 

Fexi 

6.10 

358.694 

358.6138 

60.0 ± 11.1 



a 

28 

Fexi 

6.10 

369.205 

369.1635 

60.0 ± 13.2 




29 

Fexn 

6.20 

186.867 

186.8540 

186.8870 


498 ± 64.0 


c 

30 

Fexn 

6.20 

192.390 

192.3940 


206 ±25.1 



31 

Fexn 

6.20 

193.505 

193.5090 


487 ± 55.4 



32 

Fexn 

6.20 

195.117 

195.1190 


869 ± 100 



33 

Fexn 

6.20 

196.642 

196.6400 


120 ± 18.0 



34 

Fexn 

6.20 

211.736 

211.7320 


94.2 ±29.1 



35 

Fexn 

6.20 

217.277 

217.2760 


97.7 ± 35.3 


b, c 

36 

Fexn 

6.20 

219.449 

219.4370 


73.5 ± 22.5 



37 

Fexn 

6.20 

338.273 

338.2630 

92.6 ± 11.3 


57.6 ± 9.50 


38 

Fexn 

6.20 

346.849 

346.8520 

115 ± 14.0 


77.9 ± 12.8 


39 

Fexn 

6.20 

352.112 

352.1060 

190 ±22.1 


171 ± 27.0 


40 

Fexn 

6.20 

364.494 

364.4670 

257 ± 29.5 




41 

Fe xiii 

6.25 

196.519 

196.5253 


135 ± 33.4 



42 

Fe xm 

6.25 

200.017 

200.0216 


304 ± 36.7 



43 

Fe xm 

6.25 

201.118 

201.1259 


470 ± 60.6 



44 

Fe xm 

6.25 

202.042 

202.0443 


1210 ± 138 



45 

Fe xm 

6.25 

203.164 

203.1653 


154 ± 24.2 



46 

Fe xm 

6.25 

203.820 

203.7722 

203.7957 


1930 ± 238 



47 

Fe xm 

6.25 

204.255 

204.2626 


195 ± 29.7 



48 

Fe xm 

6.25 

204.950 

204.9422 


267 ± 44.5 


c 

49 

Fe xm 

6.25 

209.623 

209.6199 


210 ±33.1 



50 

Fe xm 

6.25 

209.908 

209.9167 


225 ± 43.6 


c 

51 

Fe xm 

6.25 

213.764 

213.7686 


91.7 ±20.6 



52 

Fe xm 

6.25 

221.813 

221.8281 


218 ± 32.5 



53 

Fe xm 

6.25 

240.687 

240.6964 


160 ±66.2 


a, b 

54 

Fe xm 

6.25 

246.187 

246.2095 


246 ± 47.4 



55 

Fe xm 

6.25 

251.939 

251.9529 


332 ± 65 



56 

Fe xm 

6.25 

311.583 

311.5475 



16.4 ± 7.40 

b 

57 

Fe xm 

6.25 

312.172 

312.1748 

105 ± 16.5 

172 ± 33.8 

61.6 ± 12.0 

c* 

58 

Fe xm 

6.25 

312.897 

312.8683 

36.8 ±9.13 


33.5 ± 7.90 


59 

Fe xm 

6.25 

318.125 

318.1302 

85.9 ± 12.6 

112 ±42.8 

64.7 ± 11.6 

b* 

60 

Fe xm 

6.25 

320.807 

320.8010 

170 ± 19.8 

195 ±57.8 

118 ± 19.0 


61 

Fe xm 

6.25 

321.475 

321.4662 

31.1 ± 5.52 


28.7 ± 6.00 


62 

Fe xm 

6.25 

348.179 

348.1840 

210 ±24.3 


170 ± 27.0 

c* 

63 

Fe xm 

6.25 

359.665 

359.6445 

170 ± 19.7 
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Ion 

logT 

SERTS 
A. (A) 

CHIANTI 
A (A) 

SERTS-93 

SERTS-95 

SERTS-97 

Notes 

64 

Fe xm 

6.25 

359.852 

359.8399 

41.7 ±5.81 




65 

Fexiv 

6.30 

211.317 

211.3172 


2510 ±287 



66 

Fexiv 

6.30 

219.121 

219.1305 


748 ± 89.1 


d 

67 

Fexiv 

6.30 

220.090 

220.0849 


621 ± 84.2 



68 

Fexiv 

6.30 

252.168 

252.1996 


364 ± 73.3 



69 

Fexiv 

6.30 

257.370 

257.3941 


295 ± 96.3 



70 

Fexiv 

6.30 

264.768 

264.7889 


1270 ± 156 



71 

Fexiv 

6.30 

270.511 

270.5208 


674 ± 118 



72 

Fexiv 

6.30 

274.193 

274.2037 

1140 ± 167 

1360 ± 176 


c 

73 

Fexiv 

6.30 

289.140 

289.1508 

60.5 ±21.6 

213 ±46.4 


b*,c* 

74 

Fexiv 

6.30 

334.173 

334.1783 

645 ± 72.9 

923 ± 107 

554 ± 89.0 

c* 

75 

Fexiv 

6.30 

353.856 

353.8364 

291 ± 32.9 




76 

Fexv 

6.35 

243.771 

243.7940 


668 ± 102 



77 

Fexv 

6.35 

284.151 

284.1630 

5780 ±651 

13900 ± 1590 



78 

Fexv 

6.35 

302.322 

302.3340 



37.3 ± 10.3 

a 

79 

Fexv 

6.35 

304.873 

304.8940 



102 ± 18.0 


80 

Fexv 

6.35 

307.764 

307.7470 



43.1 ± 11.1 


81 

Fexv 

6.35 

312.566 

312.5590 

41.4 ±9.02 


51.7 ± 10.5 

a 

82 

Fexv 

6.35 

317.640 




3.60 ± 3.10 

b 

83 

Fexv 

6.35 

321.787 

321.7690 

16.3 ± 3.64 


23.4 ± 5.30 


84 

Fexv 

6.35 

324.975 

324.9750 



7.80 ± 3.60 

b 

85 

Fexv 

6.35 

327.032 

327.0330 

52.2 ± 6.64 


74.4 ± 12.1 


86 

Fexv 

6.35 

330.994 




6.20 ± 3.60 

b 

87 

Fexv 

6.35 

417.296 

417.2580 

241 ±29.5 




88 

Fexvi 

6.45 

251.053 

251.0630 


521 ±74.6 



89 

Fexvi 

6.45 

262.965 

262.9760 


858 ± 104 



90 

Fexvi 

6.45 

335.404 

335.4100 

4090 ± 470 

12000 ± 1350 

5020 ± 804 


91 

Fexvi 

6.45 

360.782 

360.7590 

2040 ± 232 




92 

Fe xvn 

6.70 

347.848 

347.8170 

9.17 ± 3.73 


17.2 ± 6.30 

b, c 

93 

Fe xvii 

6.70 

350.489 

350.4782 

17.0 ± 6.43 


18.9 ±5.90 

b*,c 


Notes. 

a CHIANTI indicates that these lines are blended; they are not used in our analysis. 
b These detections are <3tr; * indicates only one data set value is problematic. 

c Predicted intensity is >5 a low compared with observed; * indicates only SERTS-95 value is problematic. 
d Predicted intensity is >5cr high compared with observed. 


Our second DEM method is xrt_dem_iterative2, a program 
first developed (Weber et al. 2004) and tested (Schmelz et al. 
2009) for XRT data alone, but now applied more generally (see, 
e.g., Schmelz et al. 2010; Winebarger et al. 2011). The routine 
employs a forward-fitting approach where a DEM is guessed and 
folded through each response to generate predicted fluxes. This 
process is iterated to minimize / 2 for the predicted-to-observed 
flux ratios. The DEM function is interpolated using Nj — 1 
splines, representing the degrees of freedom for Nj observations. 
This routine uses Monte Carlo iterations to estimate errors on 
the DEM solution. For each iteration, the observed flux in each 
line/filter was varied randomly and the program was run again 
with the new values. The distribution of these variations was 
Gaussian with a centroid equal to the observed flux and a width 
equal to the uncertainty. 

The intensity for each spectral line observed by SERTS is 
proportional to x G(T) x DEM(T) x AT, where A is 
the elemental abundance, G(T) is the contribution function 
(erg cm 3 s _1 ) from CHIANTI, and DEM is the differential emis- 
sion measure (cnT 5 K ). The intensity for each SXT filter is 
proportional to Resp(T) xDEM(7)A7’, where Resp(T) is the 
response function (DN s -1 pixeR 1 per unit emission measure ) 
calculated from the instrument effective areas available in 
SolarSoft and the CHIANTI synthetic solar spectrum with the 
same abundance and ionization equilibrium assumptions used 


for the spectral lines. Note that the SERTS and SXT data do not 
have to be converted to the same units, as long as the proper 
contribution and response functions are used. 

The DEM curves in Figure 1 were constructed using only the 
SERTS iron lines (Table 1) and the SXT broadband filter data, 
which depend strongly on the strengths of the iron lines in the 
synthetic spectrum that is folded through the effective area of 
each filter. These curves were constructed with DEM_manual. 
Figure 2 shows the comparison of the DEM_manual curve in 
black and the Monte Carlos from xrt_dem_iterative2 in red. This 
plot is made with the SERTS-97 data set only, but the results 
apply more generally to SERTS-93 and SERTS-95. Both the 
agreement between the different methods and the tightness of 
the Monte Carlos are encouraging, showing that the DEM shape 
is indeed determined by the data. Integrating over the coronal 
temperature range for the curves in Figure 1, we find that the 
total column emission measure is EM = 1.2 x 10 28 cm~ 5 
for SERTS-93, EM = 2.8 x 10 28 cm" 5 for SERTS-95, and 
EM = 1.5 x 10 28 cm -5 for SERTS-97. 

Figure 3 shows results for the SERTS-95 data using 
xrt_dem_iterative2. The DEM in the top panel uses the SERTS 
iron lines only (without SXT). The Monte Carlos begin to devi- 
ate above log T — 6.6, indicating that the DEM is not well con- 
strained by the data. The bottom panel shows the improvement 
when SXT data are included in the analysis. The distribution 
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Table 2 

SERTS Non-iron Line Intensities (erg cm~ 2 s -1 sr -1 ) 



Ion 

log T 

SERTS 
A (A) 

CHIANTI 

A (A) 

SERTS-93 

SERTS-95 

SERTS-97 

Notes 

1 

Mg v 

5.50 

351.114 

351.0850 



16.7 ± 4.90 


2 

Mg vi 

5.65 

349.149 

349.1249 

349.1639 

80.3 ± 15.1 


55.9 ± 10.5 


3 

Mg vii 

5.80 

319.026 

319.0340 

51.9 ±8.2 


55.9 ± 10.5 


4 

Mg vii 

5.80 

367.698 

367.6780 

367.6880 

49.5 ±11.1 




5 

Mg viii 

5.90 

311.769 

311.7730 

39.5 ± 7.85 


20.9 ± 6.00 


6 

Mg vm 

5.90 

313.734 

313.7440 

67.8 ± 10.5 


71.5 ± 13.4 


7 

Mg vm 

5.90 

315.022 

315.0160 

232 ± 27.0 

144 ± 46.4 

183 ±30.0 


8 

Mg vm 

5.90 

317.017 

317.0280 

55.9 ± 10.2 

62.1 ± 25.7 

34.8 ± 8.40 

a* 

9 

Mg vm 

5.90 

338.997 

338.9840 

64.2 ± 8.48 


29.7 ± 5.40 


10 

Mg ix 

6.00 

368.093 

368.0713 

943 ± 110 




11 

Alx 

6.10 

332.784 

332.7900 

106 ± 12.7 

173 ±44.1 

1 12 ± 18.0 


12 

Si vm 

5.95 

214.757 

214.7590 


102 ± 37.6 


b, c 

13 

Si VIII 

5.95 

314.358 

314.3560 

56.0 ± 10.2 


36.8 ± 8.50 


14 

Si vm 

5.95 

316.208 

316.2180 

112 ± 17.4 


71.1 ± 12.8 


15 

Si vm 

5.95 

319.841 

319.8400 

148 ± 18.4 


91.2 ± 15.3 


16 

Si ix 

6.05 

258.077 

258.0823 


170 ± 49.2 


b 

17 

Si ix 

6.05 

292.756 

292.7594 


126 ± 49.3 


a, b 

18 

Si ix 

6.05 

292.858 

292.8092 

292.8546 


61.6 ± 44.2 


a, b 

19 

Si ix 

6.05 

296.108 

296.1135 

270 ±41.0 




20 

Si ix 

6.05 

341.971 

341.9511 

44.3 ± 6.25 


28.3 ± 5.20 


21 

Si ix 

6.05 

344.987 

344.9543 



16.7 ±4.10 


22 

Si ix 

6.05 

345.138 

345.1210 

124 ± 15.5 


65.7 ± 10.9 


23 

Si ix 

6.05 

349.871 

349.7917 

349.8602 

197 ± 23.2 


121 ± 20.0 


24 

Six 

6.15 

258.361 

258.3742 


379 ± 71.8 



25 

Si ix 

6.15 

261.043 

261.0567 


265 ± 66.9 


c 

26 

Si ix 

6.15 

347.401 

347.4026 

286 ± 32.6 


187 ±30.0 


27 

Si ix 

6.15 

356.052 

356.0381 

356.0496 

255 ± 30.0 




28 

Si xi 

6.20 

303.323 

303.3268 

2880 ± 326 

3360 ± 386 

2340 ± 370 

c* 

29 

Si xi 

6.20 

365.431 

365.4390 

65.6 ± 11.1 




30 

Sx 

6.15 

259.465 

259.4967 


216 ±66.9 


b 

31 

Sxi 

6.25 

188.663 

188.6753 


56.4 ± 14.0 



32 

Sxi 

6.25 

247.143 

247.1594 


151 ±53.3 


a, c 

33 

Sxi 

6.25 

281.364 

281.4021 


208 ± 79.8 


a, c 

34 

Sxi 

6.25 

285.820 

285.5875 

285.8226 


132 ± 44.3 


a 

35 

Sxii 

6.30 

212.088 

212.1205 


78.6 ± 16.7 



36 

Sxn 

6.35 

218.196 

218.2005 


156 ± 29.9 


b 

37 

Sxii 

6.35 

288.415 

288.4210 

136 ±35.9 

203 ± 54.7 



38 

Sxii 

6.35 

299.553 

299.5407 



34.7 ± 11.5 


39 

Sxiii 

6.40 

256.666 

256.6852 


802 ± 108 



40 

Sxiv 

6.50 

417.688 

417.6611 

67.1 ± 12.8 




41 

Arxi 

6.30 

188.821 

188.8060 


28.3 ± 10.7 


a 

42 

Arxiv 

6.55 

187.962 

187.9690 


26.3 ± 8.60 



43 

Arxiv 

6.55 

194.401 

194.3960 


18.5 ± 6.80 


a 

44 

Ca vii 

5.70 

342.435 

342.3953 

342.8179 



7.3 ±3.10 

a 

45 

Caxiv 

6.45 

186.616 

186.6103 


77.4 ± 15.5 



46 

Caxiv 

6.55 

193.872 

193.8661 


76.2 ± 10.7 



47 

Ca xv 

6.65 

200.978 

200.9719 


118 ±20.0 



48 

Ca xvii 

6.70 

192.883 

192.8532 


39.0 ± 10.7 


b 

49 

Ni xvi 

6.45 

194.026 

194.0460 


19.6 ± 3.50 



50 

Ni xvi 

6.45 

239.504 

239.5080 


156 ±61.0 


a 

51 

Ni xvn 

6.50 

249.168 

249.1856 


587 ± 81.8 



52 

Ni xviii 

6.55 

291.983 

291.9840 

123 ± 28.7 

481 ± 68.7 



53 

Ni xviii 

6.55 

320.559 

320.5660 

46.6 ±6.59 


85.3 ± 14.3 


54 

Cr xiii 

6.20 

328.257 

328.2680 

57.7 ±7.32 


62.9 ± 10.3 


55 

Crxiv 

6.25 

389.903 

389.8640 

55.9 ± 12.3 




56 

Mnxv 

6.35 

361.007 

361.0120 

29.1 ±5.56 





Notes. 

a These detections are <3<x; * indicates only one data set value is problematic. 
b CHIANTI indicates that these lines are blended; they are not used in our analysis. 

c Predicted intensity is >5cr low compared with observed; * indicates only SERTS-95 value is problematic. 
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Figure 1. DEM distributions for the active regions observed during the three SERTS rocket flights constructed with the SERTS iron lines and constrained at the 
high-temperature end by the SXT broadband data which are dominated by iron lines at these temperatures. These distributions conform to the definition of DEM as 
n^dl/dT and use the ionization fractions of Bryans et al. (2009) and the recommended coronal abundances of Schmelz et al. (2012). 

(A color version of this figure is available in the online journal.) 



og(T[K]) 


Figure 2. Comparison of the DEM_manual curve in black and the Monte Carlos from xrt_dem_iterative2 in red for the SERTS-97 data. The agreement between the 
different methods and the tightness of the Monte Carlos are encouraging, showing that the DEM shape is indeed determined by the data. 

(A color version of this figure is available in the online journal.) 


of Monte Carlos is tight, even at high temperatures, indicating 
that the DEM is now well constrained. In the following analysis, 
we use the DEM results obtained with the combination of the 
SERTS iron lines and the SXT data. 

3.2. Deriving the Electron Density from DEM Analysis 

The density is a key parameter of the coronal plasma, essential 
for trying to understand the fundamental physics of complex 
phenomena such as coronal heating and loop dynamics. Plasma 
densities are usually determined from the ratio of intensities of 
two spectral lines, ideally from the same element and ionization 
state to avoid having to know the elemental abundances and the 
ionization fractions. For example, Brosius et al. (1996) used this 


line-ratio method to calculate the electron density for the core 
of AR 7563 from several pairs of iron lines observed during 
SERTS-93. Their results showed a wide range of densities 
(1.5 x 10 9 to 2.0 x 10 10 cm -3 ), even within the same ion 
species (Fexm) where one might expect more consistency in 
the absence of systematic uncertainties in the atomic data. 

The plots in Figure 4 represent the SERTS-93 data set only, but 
the results apply more generally to SERTS-95 and SERTS-97 
as well. These panels show the predicted-to-observed intensity 
ratios for each of the iron lines listed in Table 1, for three 
different adopted values of the plasma electron density. The 
black stars represent spectral lines that are not sensitive to 
density changes (in this range), the blue squares represent 
lines whose predicted-to-observed intensity ratio decreases with 
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Figure 3. DEM results for the SERTS-95 data using xrt_dem_iterative2 with the SERTS iron lines only (top panel) and the SERTS iron lines plus SXT data (bottom 
panel). The minimum y 2 DEM is shown in black and the Monte Carlos are in red. The Monte Carlos begin to deviate above log T = 6.6 in the top panel but are tight 
in the bottom panel, which illustrates the need for the SXT high-temperature constraint (see the text). 

(A color version of this figure is available in the online journal.) 


increasing density, and red diamonds represent lines whose 
predicted-to-observed intensity ratio increases with increasing 
density. In the top plot ( n e = 2 x 10 8 cm -3 ; / 2 = 15.9), the 
squares are too high and the diamonds are too low; the reverse is 
true for the bottom plot (n e — 2 x 10 10 cm” 3 ; x 2 = 11.1). The 
middle plot (n e = 2 x 10 9 cnT 3 ; y 2 = 2.7) gives the best match 
of predicted-to-observed intensities for the three densities tested. 
The best-fit densities using this method are n e — 2.2 x 10 9 cm” 3 
with x 2 = 2.3 for SERTS-93, n e = 2.9 x 10 9 cm' 3 with x 2 = 
7.4 for SERTS-95, and n e = 2.3 x 10 9 cm" 3 with x 2 = 1.5 for 
SERTS-97. 

Using these densities, the total column emission measure 
calculated in the previous subsection, and assuming a path length 
of(l)~2 x 10 9 cm, a typical loop width, or (2) ~2 x 10 10 cm, 
a typical loop height (Bray et al. 1991), we can calculate the 
geometric density, n g , and the resulting filling factor for each 
SERTS observation: (1) n g = 2.4 x 10 9 cm” 3 and ff = 1.1 for 
SERTS-93, n g = 3.7 x 10 9 cm' 3 and ff = 1.3 for SERTS-95, 
and n g = 2.8 x 10 9 cm’ 3 and ff = 1.2 for SERTS-97. 
(2) n g = 7.8 x 10 s cm” 3 and ff = 0.35 for SERTS-93, 


n g = 1.2 x 10 9 cm” 3 and ff = 0.41 for SERTS-95, and 
n g = 8.7 x 10 9 cm” 3 and ff = 0.38 for SERTS-97. 

This illustrates the problem with the geometric method: it 
depends strongly on the estimate of the path length, a parameter 
which is not easy to pin down. For typical values, the filling 
factor can range from ~30% to 40% up to one, within the 
uncertainties. These values are consistent with the range of 
filling factors for active regions present in the literature. 

Our approach for constraining the density seemed promising 
for SERTS-93 (Figure 4) and SERTS-97, but the agreement 
between the observed and predicted values for some iron lines in 
the SERTS-95 data set was so bad that we would never get close 
to an acceptable value for x 2 if we used all the available lines. 
A number of line predictions might well be inaccurate because 
of problems with the atomic data or ionization fraction and 
there might be systematic errors for observed intensities with 
poor signal to noise. Therefore, we investigated the following 
possibility: could we get an acceptable fit for a single density 
by successive elimination of lines with the greatest predicted- 
to-observed discrepancies from the derived DEM? If we could 
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Figure 4. Predicted-to-observed intensity ratios for the SERTS-93 iron lines listed in Table 1. The points are plotted in temperature order for different values of the 
electron density: (a) n e = 2 x 10 8 cm -3 and y 2 = 15.9, (b) n e = 2 x 10 9 cm -3 and y 2 = 2.7, and (c) n e = 2 x 10 10 cm -3 and y 2 = 1 1.8. The black stars represent 
spectral lines that are not sensitive to density changes (in this range), the blue squares represent lines whose predicted-to-observed intensity ratio goes down as density 
increases, and red diamonds represent lines whose predicted-to-observed intensity ratio goes up as density increases. 

(A color version of this figure is available in the online journal.) 


achieve / 2 ~ 1 while keeping adequate temperature coverage 
with a sufficient number of density-sensitive lines, that would 
allow us to estimate la uncertainties on the density by finding 
the density values for which the total / 2 increased by unity from 
the best-fit value (i.e.. A/ 2 = Xmin + L see, e -g-> Bevington & 
Robinson 1992, pp 144-147). 

The results of this approach for the SERTS-95 data are shown 
in Figure 5. The dotted curve shows the initial x 2 distribution 
where we used all available iron lines, the dot-dashed curve 
shows the changes if we remove lines for which the predicted- 
to-observed intensity ratio was >5a from unity, the dashed 
curve shows the changes if we remove lines for which the 
predicted-to-observed intensity ratio was >3 a from unity, and 
the solid curve shows the final distribution where we remove 
even those lines with predicted-to-observed intensity ratios that 
were >2 a from unity. This method left adequate coverage for 
fitting density over the temperature range 6.0 ^ log T ^ 6.4. 
The best-fit density is given by the lowest point in the solid 
curve. The density for the SERTS-95 data set obtained using 
this method is n e = (6.9+ 2 {' 2 ) x 10 9 cm~ 3 . As noted above, the 


quoted uncertainties are la and correspond to an increase of one 
in Ax 2 . 

There is no a priori reason to think that all the EUV /X-ray 
emission from the core of a given active region would come 
from loops with a single mean electron density across a broad 
temperature regime. Rather, these DEM-weighted values are 
exactly those densities responsible for the bulk of the observed 
emission. Because the method uses a large number of spectral 
lines simultaneously, it is not heavily biased by potential atomic 
data uncertainties in any given line ratio. We suggest that 
determining the DEM-weighted density might be a powerful 
density diagnostic tool, complementary to the standard line- 
ratio method. In addition, DEM forward-fitting may provide a 
straightforward method to estimate realistic uncertainties. 

3.3. Problematic Lines and Uncertainties 

Figure 6 shows the predicted-to-observed intensity ratios for 
the SERTS-93 lines. The top panel shows the iron lines, the 
middle panel shows the non-iron lines with Schmelz et al. (20 1 2) 
abundances, and the bottom panel shows the non-iron lines with 
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Figure 5. Plots of / 2 as a function of electron density for SERTS-95. The dotted curve shows the initial y 2 distribution, the dot-dashed curve shows the changes if 
we remove lines where the predicted-to-observed intensity ratio was >5 a from unity, the dashed curve shows the changes if we remove lines where the predicted-to- 
observed intensity ratio was >3a from unity, and the solid curve shows the final distribution where we remove lines where the predicted-to-observed intensity ratio 
was >2 a from unity. The minimum in the solid curve represents the best-fit density, where the y 2 has dropped to ~1. 

(A color version of this figure is available in the online journal.) 


the abundances of the various elements adjusted to minimize 
the x 2 . The peak formation temperatures of the lines generally 
increase from left to right across the plot. Figures 7 and 8 show 
similar results for the SERTS-95 and SERTS-97 lines. 

There are many SERTS detections that are <3 a in signif- 
icance. These are labeled with a “1” in the last column of 
Tables 1 and 2. Lines that are blended are labeled with a “2” and 
are not used in the DEM analysis. Lines where the predicted-to- 
observed intensity ratios are >5cr lower than unity are labeled 
with a “3” and lines where the predicted-to-observed inten- 
sity ratios are >5cr higher than unity are labeled with a “4.” 
Adjusting the DEM curve, the electron density, or the elemen- 
tal abundance did not improve the fit. The problems with these 
lines may result from unidentified blends or poor signal to noise. 
There are also two Fe xvii lines where the differences between 
the observed and predicted values for the intensity are so large 
that we suspect that these lines may be misidentified, blended, 
or require improved atomic data. These lines were also excluded 
from the DEM analysis. 

The notes column of Tables 1 and 2 shows that there is only 
one line labeled with a “4,” the Fe xiv line at219.12lA observed 
as part of SERTS-95. It is a relatively strong, well-detected 
line and appears to be properly identified. A high predicted-to- 
observed intensity ratio could result from resonance scattering, 
but this seems unlikely (see, e.g., Brickhouse & Schmelz 2006). 
A more plausible explanation may be related to the density 
sensitivity since it does not appear to agree with the results of 
the other SERTS-95 Fexiv lines. Many more lines, however, 
are labeled with a “3” in the notes column indicating that the 
predicted-to-observed intensity ratios are too low. This implies 
that the problem may be either missing lines (unknown blends) 
or missing flux from processes that are not yet accounted for. 
This problem is more prevalent for the SERTS-95 lines, perhaps 
because this wavelength range is not as well studied. 


3.4. Relative Elemental Abundances 

We were interested in seeing how well the three DEM 
curves in Figure 1, created with iron lines alone, predicted 
the intensities for the non-iron lines. Initially, the unadjusted 
abundance values (relative to the iron abundance) from Schmelz 
et al. (2012) were used for each separate SERTS data set. The 
results of the analysis are shown in the three middle panels of 
Figures 6-8. We discuss each element separately below. 

Magnesium. The DEM curves optimized for the iron lines 
do a good job of reproducing the magnesium line intensities. 
There is one caveat, however. Many of the lines considered 
here are from the Mgv-Mgvui ionization stages, with peak 
formation temperatures of log T — 5. 5-5. 9. We have to bootstrap 
the iron-only DEM curves to slightly lower temperatures. As 
a result, there is some uncertainty in the exact placement of 
the low-temperature end of the DEM curve, and therefore, the 
exact value of the magnesium abundance required to bring the 
observed and predicted values for the intensities into agreement. 
There is, however, Fe x-Mg ix overlap for SERTS-93, each with 
a peak formation temperature of log T — 6.0, as well as Fe vm 
and Feix in SERTS-95. In summary, the SERTS data give us 
no reason to think that the recommended abundance value for 
magnesium from Schmelz et al. (2012) needs to be adjusted. 

Aluminum. The Alx line at 332.784 A is the only aluminum 
line observed, but it is listed in all three data sets. There are two 
strong detections: 8 a in SERTS-93 and 6.2a in SERTS-97. 
The SERTS-95 detection is only 4a. The line has a peak 
formation temperature of log T — 6.10, so its predicted intensity 
is well constrained by the iron-line DEMs. The SERTS-93 
aluminum abundance needs to be reduced to 64% ± 10% of the 
Schmelz et al. (2012) recommended value in order to minimize 
X 2 , whereas the SERTS-95 value must be increased to 190% ± 
70%. The 1997 value is consistent with the recommended 
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Figure 6. Predicted-to-observed intensity ratios for the spectral lines observed by SERTS-93. Top: the iron lines only; middle: the non-iron lines with the abundance 
normalization of Schmelz et al. (2012); and bottom: the non-iron lines with the abundance values adjusted to best match the total SERTS data set. 


abundance. Strictly speaking, the SERTS results themselves are 
not mutually consistent. However, we hesitate to draw strong 
conclusions based on a single line, and do not necessarily 
consider this difference to be the result of abundance variability. 

Silicon. There were more lines of silicon than of any other 
element except iron, with 11, four, and nine usable lines in the 
SERTS-93, -95, and -97 data sets, respectively. These silicon 
lines were produced from a variety of ionization stages from 
Si viii to Si xi, with a range of peak formation temperatures of 
log T — 5. 9-6. 2. The iron-only DEM curves do a good job of 
constraining all of these except for the SERTS-93 and -97 Si vm, 
but these were both bootstrapped using cooler magnesium lines 
(see above). The SERTS-93 and -97 data indicate that the silicon 
abundance should be lowered to 85% ± 10 % and 75% ± 10 % 
of the Schmelz et al. (2012) recommended value, respectively. 
The four SERTS-95 silicon lines do not agree with the scaling or 
with each other. We suspect that the SERTS-95 Si vm 2 14.76 A 
line and the Six 261.04 A line may be blended with lines that 
are not yet part of the CHIANTI database. If we eliminate these 
lines from our analysis, our results indicate that the silicon 


recommended abundance should be lowered to ~85% of its 
current value. 

Sulfur. The elemental abundance of sulfur is notoriously 
difficult to pin down because it (along with phosphorus) is the 
so-called intermediate-FIP element. Two SERTS-95 S xi lines 
at 247.14 and 281.36 A are consistent with unidentified blends 
and inconsistent with the results of the other sulfur SERTS lines. 
If we eliminate these lines from our analysis, our results indicate 
that the recommended sulfur abundance of Schmelz et al. (2012) 
is too high, and should be lowered to ~65% of its current value. 

Argon. Argon is the only high-FIP element in the data set. 
The only detections were from SERTS-95, and two of these 
were eliminated because they were weak (<)3(T detection). The 
best-fit abundance (102 ± 33) for the single remaining line is 
consistent with the recommended value of Schmelz et al. (2012) 
but does not tightly constrain it. 

Calcium. There is one weak (<3er) Cavil line in the SERTS- 
97 data set and four lines formed at higher temperatures in the 
SERTS-95 data set. The Caxiv line is not in the CHIANTI 
database and the Caxvn is blended. The remaining Caxiv 
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and xv lines are 5cr detections, but Figure 7 shows that the 
predicted-to-observed intensity ratio for the first is high and the 
second is low. The available SERTS calcium data do not appear 
to present a viable case for changing the recommended value of 
Schmelz et al. (2012). 

Nickel. The nickel lines from SERTS-93 and SERTS-95 
are consistent with the Schmelz et al. (2012) recommended 
abundance value, with best-fit results of 106% ± 10% and 
87% ± 16%, respectively. 

Chromium. There are two chromium lines in the SERTS-93 
data set and one in the SERTS-97 data set. All are consistent 
with the Schmelz et al. (2012) recommended value. 

Manganese. There is one Mnxv line in the SERTS-93 data 
set. Schmelz et al. (2012) had no results for manganese itself, but 
if we apply their algorithm for low-FIP elements to manganese, 
we find that the resulting abundance value is consistent with the 
SERTS data (see Figure 6). 


4. DISCUSSION 

4.1. Comparison of Forward-fitting with Inverse DEM Curves 

The differences between our SERTS-93 DEM curve in 
Figure 1 and the distribution published originally by Brosius 
et al. (1996, see their Figure 8(a)) were discussed in detail by 
Schmelz et al. (1999) and are summarized here. Both analyses 
used the same SERTS data, but the normalizations and the 
shapes of the two resulting curves are distinctly different. In 
our analysis, the normalization of the curve conforms to the 
definition of DEM as n 2 e dl /dT ; in addition, we have used the 
recommended coronal abundances of Schmelz et al. (2012) 
instead of the old values of Feldman ( 1 992). The main difference 
in the shapes of the two curves can be attributed to our use of 
cotemporal data from three SXT broadband filters (thin Al, 
AlMgMn, and thick Al) as a high-temperature constraint (see 
Schmelz et al. 1999). The Fe xv and Fe xvi SERTS lines as well 
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as the Ni xviii lines were in agreement with the SXT data. The 
predicted-to-observed ratios of the two weak Fexvn lines at 
347.8 and 350.5 A, however, were still too low by more than a 
factor of two and would have required an unrealistically large 
cross-calibration correction factor to bring them into agreement. 
Since this discrepancy could not be attributed to statistical or 
systematic uncertainties, Schmelz et al. (1999) suggested that 
the atomic physics for these lines be reexamined. 

There is no DEM curve for the SERTS-95 data published 
previously by the SERTS team or — to the best of our knowl- 
edge — by anyone else, to compare with our DEM curve. Our 
curve, shown by the dashed line in Figure 1, is well covered 
from log T — 5. 9-6.4, by multiple lines from each ionization 
stage of iron from Feix through Fexvi inclusive. There are 
also two Fevm lines, but one is blended and one is a <3er 
detection, so extending the DEM curve to lower temperatures 
is problematic. As in the SERTS-93 case, plasma at higher 


temperatures is nicely constrained by the SXT broadband filter 
data (see Figure 3). 

For the SERTS-97 data, after we adjusted the normalization 
of the DEM curve obtained by Brosius et al. (2000) as discussed 
above for the SERTS-93 data, our respective distributions are 
essentially identical for temperatures above log T— 5.8 (cf. their 
Figure 4 and the dash-dotted line in our Figure 1). We attribute 
this agreement to the fact that Brosius et al. chose not to use the 
discrepantFe xvii lines at 347.8 and 350.5 A mentioned above in 
their SERTS-97 DEM analysis. This agreement is encouraging 
since their method uses an automated inversion and ours uses 
two different forward-fitting approaches, and there have not 
been many direct comparisons of DEM curves from the two 
methods for real data published in the literature. 

Figure 1 shows that different active regions have different 
DEM shapes (in addition, see Warren et al. 2012; Schmelz 
& Pathak 2012). The agreement between the results from 
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DEM_manual and xrt_dem_iterative2 as well as the older 
inversion results (where available) described above gives us 
confidence that the SERTS plus SXT data together provide a 
robust determination of the curve, and that the differences are 
not attributable to possible idiosyncrasies in the DEM methods. 
Some of these differences are most likely related to the core 
heating of the region, whether it is steady (Winebarger et al. 
2011) or nanoflare (Tripathi et al. 2011) heating, for example. 
With these new DEM techniques and the improved atomic 
data available in CHIANTI v.7, we expect to make significant 
progress in obtaining the more accurate physical quantities 
needed to understand such important physical phenomena as 
active region heating. 

4.2. Absolute Elemental Abundances 

From a survey of absolute coronal abundance values pub- 
lished in the literature, Schmelz et al. (2012) found that the 
abundances of elements with low FIP were enhanced over 
the photospheric values reported by Caffau et al. (2011) and 
Lodders et al. (2009) by a mean value of 2.14 ± 0.46 (rather 
than the factor of four proposed by Feldman 1992); and the 
abundances of high-FIP elements were depleted by a mean fac- 
tor of 0.71 ± 0.14 (rather than 0.25 as proposed by Meyer 
1985). The resulting set of recommended coronal abundances 
of Schmelz et al. (2012) is approximately the geometric mean 
of the Feldman (1992) and the Meyer (1985) abundance sets. 

Absolute abundances in the corona are still a subject of some 
controversy. For example, the measured differences between 
Raymond et al. (1997) and Feldman et al. (1998) for streamer 
abundances may be the result of gravitational settling. The 
different iron abundance values obtained by White et al. (2000) 
using Solar and Heliospheric Observatory CDS iron lines and 
thermal bremsstrahlung emission from the hydrogen continuum 
measured by the Very Large Array (VLA) are a factor of two 
higher than those of Zhang et al. (2001) who did a similar 
analysis using EIT and the VLA. These results relied heavily 
on the instrument cross-calibrations, but there is no way to tell 
which was right (or even which was better.) The one thing 
we seem to know for certain is that coronal abundances differ 
from their photospheric counterparts in a way that depends 
primarily on FIP, although under certain conditions, additional 
physical mechanisms may need to be taken into account (see, 
e.g., Schmelz 1993). 

There is also evidence for abundance variability in flares 
(e.g., Sylwester et al. 1984), from active region to active 
region (e.g., Strong et al. 1988), and between and within 
different kinds of coronal structures (e.g., Widing & Feldman 
1989). Some of the apparent variations reported in these older 
papers might be attributable to instrument calibration problems, 
atomic physics uncertainties, or temperature effects. A deeper 
understanding of these possible variations might result from our 
improved understanding of the thermal structure of the corona, 
better multithermal analysis techniques (e.g., DEM_manual and 
xrt_dem_iterative2, with high- and low-temperature constraints 
as appropriate), and most importantly, updated atomic modeling 
for coronal emission lines now available in CHIANTI v.7. In 
fact, we tried to do this very analysis of SERTS plus SXT data 
about 10 years ago with atomic data from CHIANTI v.2 and 
the cruder DEM techniques we had available at the time. The 
results were sufficiently disappointing that we put the project 
on hold until about six months ago. 

Although it is impossible to provide a single model satisfying 
all observations, and although it is vital to account for the 


possibility of abundance variability, it is often useful to begin 
a new research project (like this one) with a given set of 
abundances. The results presented here support use of the 
recommended coronal abundances of Schmelz et al. (2012), 
which were derived from work by several groups using various 
types of data from different instruments on numerous spacecraft, 
as the default values in future investigations. 

4.3. FIP Effect 

The SERTS results reported herein suggest that two small 
changes might need to be made to the Schmelz et al. (2012) 
coronal abundance values. Silicon could be lowered to ~85% 
of its current value, which is within the 1 a uncertainty. Sulfur, 
on the other hand, might need to be lowered to ~65% of its 
current value. The lower sulfur abundance is consistent with 
the spectroscopic measurements from Veck & Parkinson (1981) 
and Fludra & Schmelz (1999), but at odds with the SEP data 
from Reames (1995) and the solar wind data from von Steiger 
et al. (2000, 2010). Abundance results for this intermediate-FIP 
element are difficult to pin down, and further analysis will be 
required to sort out this problem. 

One thing the SERTS data do not provide, however, is strong 
observational evidence of the traditional FIP step, i.e., that the 
coronal-to-photospheric abundance ratios of low-FIP elements 
seem to be enhanced relative to those of high-FIP elements. 
This is because the only high-FIP coronal detections in the 
data sets were of the argon lines from SERTS-95, which were 
simply too weak to provide strong constraints on the coronal 
abundance. The other SERTS high-FIP detections were for 
transition region lines of neon and oxygen. These lines were 
too cool to be included in our DEM analysis, which was 
constrained by the SERTS iron lines (Table 1) and the SXT 
data. Strong coronal lines from high-FIP elements tend to have 
higher energies than the SERTS lines, so future joint EUV-X- 
ray spectroscopy observations would be required to address this 
issue satisfactorily. 

5. CONCLUSIONS 

We have used the set of recommended coronal abundances 
from Schmelz et al. (2012) for our multithermal analysis of 
active region data. We feel that these abundances provide the 
best default characterization of the low corona when specific 
information on the composition of the coronal plasma under 
study is lacking. It gives the mean values of a large number of 
spectroscopic, particle, and solar wind observations from many 
different sources, and it typically reduces the uncertainties by a 
factor of two from either extreme of the two standard coronal 
abundance sets (Meyer 1985; Feldman 1992). The absolute 
normalization of the abundances is critical for determining 
the amount of emitting plasma, the radiative losses, and all of 
the physics depending on those quantities (such as the heating 
requirements and the energy budget). 

We have used two iterative forward-fitting techniques to 
determine the DEM for the core structures in three active 
regions, using high-resolution EUV spectra from three SERTS 
rocket flights, in conjunction with high-temperature constraints 
from cotemporal broadband soft X-ray data from Yohkoh SXT. 
Elemental abundance uncertainties were not a major factor in 
deriving the shape of the distribution because we chose to work 
only with iron lines in the SERTS data, and iron lines dominate 
the SXT filter responses at active region temperatures. 
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For the radiometrically calibrated SERTS-97 spectrum, the 
DEM curve we obtained with forward fitting is virtually iden- 
tical to the curve derived from an inversion method by Brosius 
et al. (2000), for temperatures above log T = 5.8. This agree- 
ment is encouraging, especially since there have been few direct 
comparisons of DEM curves obtained from the two methods for 
real data published in the literature. 

For each set of active region spectra, we were able to establish 
a characteristic electron density of the emitting plasma in the 
temperature range 6.0 ^ log T ^ 6.4 from the DEM- weighted 
density technique. We considered electron densities ranging 
from 10 s to 10 10 cm” 3 and found that the lowest y 2 values 
were achieved with n e = 2.2 x 10 9 cm” 3 with y 2 = 2.3 for 
SERTS-93, n e = 2.9 x 10 9 cm” 3 with y 1 = 7.4 for SERTS-95, 
and n e = 2.3 x 10 9 cm” 3 with y 1 = 1.5 for SERTS-97. In 
each case, an electron density significantly higher or lower 
than this best value made it increasingly difficult to derive an 
emission measure distribution with good agreement between 
theoretical and observed line intensities. We went one step 
further for SERTS-95, which had a number of bad results. We 
systematically eliminated these lines until we could get y 2 ~ 1, 
giving n e = (6.9* 2 5 6 2 ) x 10 9 cm” 3 . 

CHIANTI is a collaborative project involving the NRL 
(USA), the Universities of Florence (Italy) and Cambridge 
(UK), and George Mason University (USA). Solar physics 
research at the University of Memphis is supported by a Hinode 
subcontract from NASA/SAO as well as NSF ATM-0402729. 
The impetus for this work was a pilot study funded by a small 
Yohkoh Guest Investigator grant. 
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